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Abstract 



> 

We study the resonant interaction of charged particles with a gravitational 
| wave propagating in the non-empty interstellar space in the presence of a uni- 

form magnetic field. It is found that this interaction can be cast in the form 
of a parametric resonance problem which, besides the main resonance, allows 
for the existence of many secondary ones. Each of them is associated with 
a non-zero resonant width, depending on the amplitude of the wave and the 
energy density of the interstellar plasma. Numerical estimates of the particles' 



bJT), energisation and the ensuing damping of the wave are given. 

> : 

^ : 1 Introduction 

Despite the numerus efforts made up today to detect gravitational waves, there is no 
convincing evidence for their existence (Thorne 1987). This is due to the fact that not 
only their amplitude is very small (Smarr 1979), but it is highly possible that some 
kind of damping mechanism operates on them as they travel through space (Esposito 
1971, Macedo and Nelson 1983, Papadopoulos and Esposito 1985). This damping 
may originate in the interaction of the gravitational wave with the interstellar matter 
(Macedo and Nelson 1990, Varvoglis and Papadopoulos 1992). 

In a recent paper (Kleidis et al 1993, which hereafter is referred to as Paper 
I) the problem of the interaction of a charged particle with a gravitational wave, 
in the presence of a uniform magnetic field, has been modelled as a Hamiltonian 
dynamical system. The corresponding analysis was carried out for various directions 
of propagation of the wave with respect to the magnetic field. It was found that, in 
the oblique propagation, diffusive acceleration of the particle, due to secular energy 
transfer from the wave, could lead to a damping. 
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The most important results, however, came out from the parallel propagation case 
where the dynamical system is trapped at an exact resonance between the Larmor 
frequency of the particles and the frequency of the wave. In this case a phase lock 
situation appears (e.g. see Menyuk et al 1987), leading to an "infinite" acceleration 
of the particle and, consequently, to a non-trivial damping of the wave. The zero 
probability problem of the exact resonance was waived out in a more recent paper 
(Kleidis et al 1995, which hereafter is referred to as Paper II), by considering that 
the propagation of the gravitational wave takes place in a space filled with plasma, 
which results in a dispersion of the wave (Grishchuk and Polnarev 1980). In this case, 
resonant phenomena occur also in the quasiparallel case, i.e. propagation at a small 
angle with respect to the direction of the magnetic field, t? < 5°. 

In the present paper we perform an elaborated study of the resonant interaction 
between charged particles and a linear polarized gravitational wave, which propagates 
in a non-empty space, parallel (or quasiparallel) to the direction of a uniform mag- 
netic field B = B e z . The interstellar plasma is represented by a collisionless gas of 
particles, where by collisionless we simply mean that the mean time between succesive 
collisions of the particles is much larger than the period of the gravitational wave. 
As we show, in this case, the generalized coordinate x 1 obeys a Mathieu differential 
equation (Abramowitz and Stegun 1970). As long as the resonance condition given 
by Eq. (23) of Paper II is fullfiled, this equation corresponds to a dynamical system 
which is trapped at a parametric resonance, where any external action ammounts to 
a time variation of the frequency parameter (Landau and Lifshitz 1976). In this case, 
the system's equilibrium at rest (re 1 = 0) is unstable. Any deviation from this state, 
however small, is sufficient to lead to a rapidly increasing displacement x 1 . Then, 
the corresponding generalized momentum 7Ti satisfies a modified Mathieu equation, 
resulting in an exponential increase of the perpendicular energy, I±, of the particle. In 
this case, also, there exists a large number of secondary resonances of non-zero width 
(Bell 1957), at which the dynamical system may be trapped. Their exact location is 
given in terms of the energy parameter J . The fact that the total measure of the 
resonant widths is non-zero increases the effectiveness of the interaction mechanism 
under study, both in the process of the accelaration of particles and in the damping 
of the wave. 

2 Parametric resonance 

We consider the non-linear interaction between a charged particle and a linear po- 
larized gravitational wave propagating in a non-empty space. The interaction takes 
place in the presence of a uniform and static, in time, magnetic field, B = B e z , 
which does not interact with the gravitational wave. The motion of a charged par- 
ticle in curved spacetime is given, in Hamiltonian formalism (Misner et al 1973) by 
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the differential equations 



dx» _ OH dn^ OH 

~d\ ~ ' ~d\ ~ ~dx^ ( ' 

where 7r M are the generalized momenta (corresponding to the coordinates x M ) and the 
"super Hamiltonian" , H, is given by the relation 

H = \<T K - eA,) (*„ - eA u ) ee \ (2) 

(in a system of geometrical units where H — c — G — 1). In Eq. (2) g^ v denotes the 
components of the contravariant metric tensor, which are defined as 

with rj^ u = diag(l, —1,-1,-1) and \ W\ < 1. is the vector potential, correspond- 
ing to the tensor of the electromagnetic field in a curved spacetime F^. The mass of 
the particle is taken equal to I. For the specific form of the magnetic field we take 

A = A 1 = A 3 = , A 2 = -Box 1 

According to Paper II the problem depends on four parameters: The normalized, 
dimensionless amplitude of the gravitational wave, a, the angle of propagation with 
respect to the direction of the magnetic field, i?, the dimensionless frequency, v, (which 
is the ratio between the frequency of the wave ui and the Larmor frequency of the 
particles Vi) and the small parameter (3 ~ f?/^ 2 , which is a dimensionless constant 
depending on the energy density of the interstellar plasma {(3 <C 1) (e.g. see Grishchuk 
and Polnarev 1980, Kleidis et al 1995). 

Resonant phenomena become dominant, in particular, when the gravitational wave 
propagates parallel to the direction of the magnetic field, i.e. i? = 0°. In this case, 
the non-zero components of the metric tensor are 

= diag ( 1, — 7 , l — , -1 ] 

\ 1 — acosk^x^ 1 + acosk^ J 

and the super Hamiltonian, given by Eq. (2), is finally written in the form 

H = l (2 „ _!)(!_ ^ p )p 

- f(v, S)hh ~~ 2,nff* - \ 1 - clnfi (3) 
where, for the sake of convenience, we have set 

2u- 1 
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Since 8° is a cyclic coordinate, the corresponding generalized momentum Jo is a con- 
stant of the motion and, according to Paper II, I 7^ 0. Now, the equation of motion 
for 9 3 is readily solved to give 

where 7r/2 is the initial phase. Since a <C 1, in what follows we keep terms up to first 
order in a. Using the approximation 

1 =f asinO 3 



1 ± asin9 3 

the geodesic differential equation of motion for x 1 is written in the form 

— + «/(!/, P)I sin /?) J A] — + x 1 = (5) 

where we have also used the equation of motion for i\\ and the fact that 

dx 1 

hi = — + eAi 
dX 

in which, for the specific form of the magnetic field, Ai — 0. 

We consider a transformation of the independent variable, A, of the form s = s(A), 
where s(A) is chosen as to satisfy the differential equation 

d^ s <7s 

^ + af(u,(3)I sin[f(u,(3)I X] — = (6) 

The condition (6) guarantees that the resulting, in terms of s, equation of motion 
for x 1 will continue to represent the differential equation of a geodesic (Papapetrou 
1974). Eq. (6) may be solved in the first order to a, to give 

S (A) ~ n (a + j^fMf^ PW]^ (7) 

where the Larmor frequency, Q, appears as an integration constant. Now, in terms 
of s(A), Eq. (5) may be written in the form 

(fx 1 

— — + Q 2 ( y l-2acos[f(u,p)I ns})x 1 = (8) 
ds 2 

where we have also used the fact that This is a Mathieu equation (Abramowitz 

and Stegun 1970) which is closely related to the problem of parametric resonance in 
dynamical systems (Landau and Lifshitz 1976). According to it, the external action 
(in our case the gravitational wave) amounts only to a periodic time variation of the 
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frequency parameter of the unperturbed system (the Larmor frequency, Q). Indeed, 
in this case, the generalized frequency 

a 2 (s) = n 2 [l-2a cos( 7 s)] (9) 

differs only slightly from the constant Q and is a simple periodic function. In Eq. 
(9) 7 = f(v,{3)I Q denotes the frequency of this periodic function. It can be proved 
(Landau and Lifshitz 1976) that the main parametric resonance in a dynamical system 
appears when 7 = 2Q. 

It has been shown (Kleidis et al 1995) that the system gravitational wave + 
charged particle is trapped in a resonance when the resonant condition, given by Eq. 
(23) of Paper II, is satisfied, that is when 

/?)/o = 2 (10) 

which gives 7 = 2VL. Therefore, once the condition (10) is satisfied, our dynamical 
system is actually trapped in a main parametric resonance. In general, there is a 
region of instability of non-zero measure, the so called resonant width (Landau and 
Lifshitz 1976), around the exact value 7 = 2f2, inside which the parametric resonance 
mechanism operates, which is 

2^-^Ae < 7 < 2Q + -Ae (11) 

The range (Bell 1957) of the resonant width is given, in our case, by Ae = 2afl. 

Apart from the case where 7 is close to 2fl, parametric resonance also occurs when 
the frequency of the parameter a(s) is close to any value 2Q/n, where n is a natural 
number (Landau and Lifshitz 1976). Therefore, except from the primary resonance, 
there is also a large number of secondary ones at which the dynamical system under 
consideration may be trapped. Each of them also has a non-zero resonant width, the 
exact form of which, in terms of is given by (Bell 1957) 

V a [(n-1)!] 2 ^ [iZ) 

that is, it decreases rapidly with increasing n, as a™. 

In our case it is more convenient to express both the possition of each of the 
secondary parametric resonances and the corresponding resonant widths in terms of 
the constant of the motion J , which corresponds to a measure of the total energy of 
the dynamical system. Since, according to Paper II, the exact resonance is achieved 
for 

/ - I* - 2 

the dynamical system will be trapped at a secondary one, of order n, when 

h = 11 (13) 

n 
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For each of these secondary resonances the corresponding resonant width is given by 

A/ (n) = 1 

° f(u,p) n 

or else 

where / = f{v,(3). Using Stirling's formula, Eq. (14) is simplified to 



If, furthermore, we set 



e2 

X = j« < 1 



the resonant width of a parametric resonance of order n is finally written in the form 



A/^ = i- *2 (16) 
7r/ rr 



In this case, we may, also, calculate the total measure of the resonant widths, in terms 
of I , which corresponds to the probability of the dynamical system to be trapped in 
a resonance. The total measure of the resonant widths is written in the form 

^ A/ °^„?> (17) 

As long as x < 1 ; the series on the r.h.s of Eq. (17) satisfies D' Alembert's crite- 
rion (e.g. see Gradshteyn and Ryzhik 1965) and therefore converges absolutely and 
uniformly. In this case the total measure of the resonant widths can be written as 

E A/ o = 4^3^(1, 1, 1, ;2, 2; X ) (18) 

where we have taken into account the properties and the series expansion of the 
generalized hypergeometric series 3-^2(1, 1, 1; 2, 2; x) ( e -g- see Erdelyi et al 1953). If 
we keep only terms of first order with respect to a and f3, Eq. (18) is written in the 
form 

e 2 

VA/ = — a + 0(a[3) + 0(a 2 ) (19) 
where e is the basis of the natural logarithms. 
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3 Numerical results 



The most interesting case of interaction between a gravitational wave propagating 
in a non-empty space and a charged particle is the quasiparallel case, in which the 
propagation takes place at a small angle with respect to the direction of the magnetic 
field. 

Since, in this case, the phase velocity of the gravitational wave is greater than 
the velocity of light (Grishchuk and Polnarev 1980), it is probable that its projection 
along the z-axis, v z = v p hCos$, to be exactly equal to c. This leads (see Paper II) to a 
combined interaction (both chaotic and resonant) between the particle and the wave, 
resulting in an overall decrease in the parallel momentum J3 (which corresponds to 
an overall increase of the perpendicular one Ii) as a function of s (see Fig. la). 

Fig. la and b: Plots of ^3 versus s (the affine parameter) for a numerically integrated 
trajectory with $ = 0.5°, v = 52.01, (3 = 10~ 5 and a = 0.002. Notice the overall decrease 
in ^3 when (a) Iq = Iq = 0.038, in contrast to the purely chaotic behaviour when (b) 
Iq = 0.0385 which lies outside the corresponding theoretical resonant width AIq = 7 ■ 10 -5 . 

Since our problem corresponds to a parametric resonance, except from the primary 
resonant condition, for Iq = Iq, there is also a large number of values of the constant of 
motion I , -^0 = Io/ n ( n — 1? 2, 3, ...), leading to secondary resonances, each of which 
has a non-zero resonant width. One expects that the above mentioned theoretical 
results in the quasiparallel propagation case should continue to hold around each of 
these secondary resonances, provided that the values of Iq lie within the corresponding 
resonant width, AI^ (in connection to this see Fig. lb). 

To check for the validity of these theoretical estimates, we integrate numerically 
the equations of motion of a charged particle in a gravitational wave for several values 
of the constant of motion I , when Iq = 0.038. Some of the corresponding results are 
presented in Figs. 2 and 3. According to these, when Iq = Iq/u, we do observe the 
expected decrease in I 3 , which corresponds to the exponentially increasing brantch of 
ii, and implies that our dynamical system is trapped in a parametric resonance. 

Fig. 2a and b: Plots of I3 versus s for a numerically integrated trajectory with $ = 0.5°, 
v = 52.01, = 10~ 5 and a = 0.002 when (a) I = iJ/2 and (b) J = JJJ/5 

The fact that the interaction of a charged particle with a gravitational wave in the 
parallel and/or the quasiparallel case corresponds to a parametric resonance, may 
impose noteworthy implications to the problem of the interaction between charged 
particles and a gravitational wave, increasing both the efficiency and the probability 
to achieve this particular interaction mechanism. 

Fig. 3a and b: Similar to Figs. 2, except that (a) I = Ifi/9 and (b) Iq = Iq/12 

As a gravitational wave propagates in the non-empty interstellar or intergalactic 
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space, it encounters several clouds with various characteristic temperatures and den- 
sities, while at the same time, it intersects the cosmic magnetic fields at various angles 
(Hillas 1984). In each case corresponds a different set of values of f3, v and °&. One 
expects that, at least for some of these sets, one of the resonant conditions is satisfied 
(since the corresponding resonant widths are non-zero), leading to a resonant inter- 
action between a charged particle and the gravitational wave. According to Paper II, 
the result of this interaction is the acceleration of the particle at very high energies. 
Since the total energy Iq is constant, the dynamical system (gravitational wave and 
charged particle) is isolated. Therefore any energisation of the particle corresponds 
to a damping of the wave. 

It is worth to note that the decrease of 7 3 in Figs. 2 and 3 is more pronounced for 
large values of n. This is in contrast to what one may have expected, since the higher 
order resonances are considered to be weaker than the lower order ones. However this 
result is only phenomenological and of no physical importance, since it has to do only 
with the slope of the straight line I3 — I±. In complete correspondence with Paper II 
we may find, from the condition H ~ 1/2, that, in momentum space, the particles 
move along the straight line 



2 V/o 



2u-l, 1+ I f(u 7 P) 1 fl 



v 2 y r ' 2 4 Wo 



(20) 



and since, at every n* h -order parametric resonance we have I = Iq/u, Eq. (20) can 
be written in the form 



71 T 

h = ~~ 2 h + 



12 "-W-? 



(21) 



_n v 2 4 

We note that for n = 1 (the primary resonance) Eq. (21) is reduced to the corre- 
sponding one of Paper II [Eq. (28)]. In the present case, as n — > 00, the straight 
line h(h) tends to become parallel to the J 3 -axis. This simply means that a small 
variation in Ii leads to a large variation of J 3 , which is exactly the case shown in Figs. 
2 and 3. 

However, there is another numerical result which may be of physical importance. It 
is related to the fact that, for several of the secondary resonances, resonant phenomena 
are observed even for values of I lying outside of the corresponding, theoretically 
estimated, resonant widths [Eq. (14)]. It should be noted that this purely numerical 
result is not observed in the case of the primary resonance, I = Iq. Some of these 
results are shown in Figs. 4. 

In Fig. 4a we present a plot of I3 versus s for a numerically integrated trajectory 
with I = f + 10~ 4 . We see that an overall decrease in J 3 occurs although, in this 

case, A4 2) ~ 10~ 7 . The same is also true for the case presented in Fig. 4b, where 
I = + 10~ 6 , far away from the corresponding resonant width, which, in this case 

is A/q 9 -* ~ 10~ 27 . Therefore, it seems that the numerically estimated resonant width, 
SIq^ , around each of the higher-order parametric resonances are many orders of mag- 
nitude larger than the corresponding theoretical one. 
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Fig. 4a and b: Plots of I3 versus s for a numerically integrated trajectory corre- 
sponding to a problem of parametric resonance which is (a) similar to the one of Fig. 
2a, but with Iq = 0.0191 = -f" + SIq and (b) similar to the one of Fig. 3a, but with 
Jo = 0.004221 = §- + 51^. In any of these cases SI^ > A/J n) , but also 51^ <C £ AJ . 



This behaviour results probably from the fact that, because of the small value of 
Iq, the secondary resonances are closely spaced. Therefore, even a small diffusion 
process due to the chaotic motion, which is always present in this case, could proba- 
bly bring the value of Jo into a resonance, even if, initially, it was outside of it. Indeed, 
as we may see from Figs. 4, the resonant behaviour becomes evident only after a con- 
siderable time interval. This time interval should correspond to the correlation time 
of the diffusive process (Farina et al 1993). Therefore, our system may be a lot more 
efficient, with respect to the resonant interaction procedure, than the corresponding 
theoretical one. Notice also, that these phenomena are not present in the case of the 
primary resonance, since the corresponding value of I lies far away from every of the 
secondary ones and hence, it possesses a well defined resonant width. For this reason, 
in this case, theoretical and numerical results are similar. 



4 Exact solutions 

In Paper II it has been shown that, when the gravitational wave propagates exactly 
parallel to the direction of the magnetic field (the z-axis), we may solve the equa- 
tions of motion of the charged particle. Their solution reveals a phase-lock situation 
(Menyuk et al 1987, Karimabadi et al 1990) resulting in an "infinite" acceleration 
of the particle along the x-axis. However, there, the solution was obtained through 
an averaging technique (e.g. see Lichtenberg and Lieberman 1983). Since the prob- 
lem is identified as that of a parametric resonance, in which the equation of motion 
along the x-axis corresponds to a Mathieu equation (Gradshteyn and Ryzhik 1965, 
Abramowitz and Stegun 1970), we may find exact solutions without resorting to any 
averaging technique. 

In the parallel propagation case, 1? = 0°, the dynamical system is trapped in a 
resonance when the condition (10) is satisfied. In this case the equation of motion 
(8), for the generalized coordinate x 1 , is written in the form 

(Pr 1 

— - + tt 2 [l - 2a cos(2tts)] x 1 = (22) 
as 2 

or else 



^ + Vl 2 [l - 2a cos(2£)] x 1 = (23) 

where we have set £ = Qs. This is a Mathieu equation with a = 1 and q = a 
(Gradshteyn and Ryzhik 1965, Abramowitz and Stegun 1970). Its solution is therefore 
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given in terms of the Mathieu functions ce\r and se±, associated with even and odd 
periodic solutions respectively. Each one of these solutions has a period of n or 2n. In 
the present paper we consider periodic solutions of period 2-rr. It can be shown that, 
for a given point (a, q) in the parameter space, there can be at most one periodic 
solution of period 2n (Abramowitz and Stegun 1970) and therefore our solution is 
unique. The general solution of Eq. (23) is written in the form 

x 1 (S)=Ace 1 (S,a) + Bse 1 (S,a) (24) 

where A, B are constants. 

For a < 1 the Mathieu functions ce\ and se\ may be expanded in power series 
of a. Considering that A = C = B, in the linear approximation to a (e.g. see 
Abramowitz and Stegun 1970, Eq. 20.2.27, p. 725), Eq. (24) becomes 

x 1 = V2Csin + - a^Csin Q + 3f ) (25) 

In this case, using the normalization conditions of the Mathieu functions (e.g. see 
Gradshteyn and Ryzhik 1965, Eqs. 6.911.3, 6.911.5. and 6.911.7), we may determine 
the exact value of the constant C. The mean square displacement of x l is defined as 

< (^) 2 >= ^ j 2 J[ xl ^] 2 ^ = C2 ( 26 ) 

so that 



C = y/< (*T > (27) 
Now, Eq. (25) is written in the form 

x\s) = V2 ^< Or 1 ) 2 >sin Q + to) - a^^/< (x 1 ) 2 >sin Q + 3to) (28) 
and therefore, for a — 0, initially (s = 0) we have 



Xl (S = 0) = y/< > (29) 

corresponding to a sort of a Brownian motion (a not unexpected result). 

However, the most important results come out from the solution of the equation 
of motion for the generalized momentum 7Ti 

*H = x l (30) 

dX 1 — asinO 3 

To solve this equation we follow the same procedure as in the case of the generalized 
coordinate x 1 . We then find that the perpendicular momentum tti obeys a modified 
Mathieu equation (e.g. see Abramowitz and Stegun 1970), of the form 



d ii\ 

ds 2 



- Q 2 [1 - 2acosh{2Qs)] m = (31) 
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Its general solution is given in terms of the radial Mathieu functions Ce\ and Se\ 
(Gradshteyn and Ryzhik 1965, Abramowitz and Stegun 1970) 



7Ti(s) = A Cei(Qs, a) + B Se^Qs, a) — A cei(ifls, a) — iB sei(ifls, a) 



which, in the first approximation to a, is written in the form 



7Ti = A 



a 



cosh(£ls) — — cosh(3Qs) 
8 



B 



a 



sinh{Vts) — —sinh(3Qs) 
8 



(32) 



Compatibility of these results with those of Paper II requires that, in this case, we 
have to set B = thus obtaining 



7Ti 



A 



a 



cosh(Qs) — —cosh(3£ls) 
8 



(33) 



From Eq. (33) it is evident that there is an exponential increase in the perpendicular 
momentum of the particles, which, in the zeroth order approximation to a, corre- 
sponds to the exponential brantch of Eq. (27) of Paper II. We also note that, in the 
exact solution case, an additional, secularly increasing term, ~ a cosh(3QX), arises. 
The corresponding perpendicular energy, ii(A), in the first order approximation to a, 
is given by 



ii ~ 2° 0s2 l^ 77 ''^) cosh(a\) 



a 



1 H — cos(3X) cosh(3a\) 
2 



+ 



(34) 



where we have used Eq. (7). Eq. (34) is in complete agreement to the results of 
Papers I and II. We see that the exponential increase in the perpendicular energy 
is, furthermore, modulated by periodic functions of time, which were smoothed out 
by the averaging technique in the previous Papers. It is important to point out 
that the phase-lock situation breaks down when l\ becomes large enough, since then 
the perturbation term of the Hamiltonian, Hi, becomes comparable to Hq and our 
approximation, in terms of the small parameter a, is not valid any more. 

Based on these results we may give an estimate of the absorption power per unit 
area of the gravitational wave, as a function of the proper distance from its source to 
Earth, for typical parameters of the interstellar gas. 

In physical units, the energy density gained by the charged particles due to their 
resonant interaction with an incident gravitational wave, within a proper time interval 
AA, is 



AE, 



gained 



Ahmpc 2 ( — J n to t 



(35) 



where mp is the proton's mass, (n a ct/ntot) is the ratio of particles which interact 
resonantly with the gravitational wave, n to t is the particles' number density of the 
interstellar matter and AIi is the energy change within the interval AA, measured in 
dimensionless units. 
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Using Eq. (29) of Paper II, also translated in physical units, for v — 1 we obtain 




AE gained = ^^hm P c\j 2nkBT P n tot (36) 

Accordingly, we consider that the interstellar space is filled with a non-relativistic 
perfect proton gas (Polnarev 1972), for which 

P = lQ n QJL]llt!L ) p = m P n tot (37) 
ur mpc z 

We take n to t = lparticle/cm 3 and T = 10K as typical mean values for the in- 
terstellar space. Then, we find that the total energy density gained by the charged 
particles due to their resonant interaction with the gravitational wave, within the 
proper time interval AA, is 

AE gained = 1.61 x 10- 39 ^ (^|) (38) 

Since, on the other hand, the system wave + particles is isolated, Eq. (38) corresponds 
also to the total energy density lost from the gravitational wave, due to the interaction 
mechanism under consideration. 

The corresponding gravitational energy flux (power per unit area) lost in the same 
proper time interval is 

AF lost = 1.21 x 1(T 29 ^ (-^-) (39) 

ur cm z sec 

Using Eq. (27) of Paper II for All, we integrate Eq. (39) over A, to obtain the total 
loss of gravitational energy flux during the whole proper time interval < A < A^, at 
which the wave travels from its source to Earth. Taking into account that for every 
realistic proper time interval Ag we have a As <C 1, we obtain 

F lost = 6.1 x l(T 30 ^a 2 A! (40) 
where A^ is normalized to ^. In physical units, Eq. (40) is written in the form 

F^ = 6.1x^a 2 (^) 2 A| (41) 

Now, A^ is measured in sec while uj and Q in Hz. However, since v — 1 we have 
uj = Q. Therefore, Eq. (41), in terms of the proper distance te between the source 
of the gravitational waves and the Earth, is finally written in the form 

F^ = 6.1xl0- 30 ^a 2 r| (42) 

where r E is measured in cm. 
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The total gravitational energy flux expected in the neighbourhood of the Earth, 
if we do not take into account any interaction with the interstellar matter (Misner et 
al 1973, Thorne 1987), is 



c 



3 



expected 



wV = 8.1xlOVa 2 (43) 



16ttG 

Dividing by parts Eqs. (42) and (43) we obtain 

= 0.753 x 10- 66 (—) 6 (—Y (44) 

Fexpected V U ) \Cm ) 

Eq. (44) gives the percentage of the flux lost due to the resonant interaction between 
a gravitational wave and a distribution of charged particles, at a proper distance te 
from the source, over the theoretically expected gravitational energy flux, at the same 
distance, in the absence of resonant damping mechanisms. 

It would be interesting to give some numerical examples of the flux damping, 
for realistic sources of gravitational waves in our galaxy. As those we choose the 
gravitational systems of bright NS/NS binaries (Thorne 1995). The gravitational 
waves generated by these sources have very low frequencies, / = 10~ 4 Hz. Since 
v — 1 we have f2 = 0.63 x 10~ 3 Hz for protons, which is true for a magnetic field of 
strength B ~ 0.5 x 10~ 6 G, a quite reasonable value for the interstellar space (Hillas 
1984). 

For gravitational waves generated from a binary NS/NS source at a distance te = 
10 Kpc from Earth (the center of our Galaxy), Eq. (44) gives 



lost 



Fexpected 



0.0116 = 1.16% 



while, for the same source at a distance of r E = 30 Kpc (the other edge of our Galaxy), 
Eq. (44) gives 

Flost = 0.1044 = 10.44% 

" expected 



5 Discussion and conclusions 

In the present paper we discuss the resonant interaction of a sinusoidal plane po- 
larized gravitational wave with the interstellar matter, in the presence of a uniform 
and static in time magnetic field. This sort of interaction arises in the parallel and 
quasiparallel directions of propagation of the gravitational wave with respect to the 
dynamic lines of the magnetic field. In any of these cases of propagation and under 
certain resonant conditions, the trigonometric term in the corresponding Hamilto- 
nian function becomes stationary. The dynamical system is therefore trapped in a 
resonance (Kleidis et al 1995). 
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Under an appropriate transformation of the affine parameter, which guarantess 
that the particles will continue to move along geodesies, we may show that this non- 
linear interaction problem corresponds to a parametric resonance. In this case the 
external action (the gravitational wave) amounts to a time variation of the frequency 
parameter of the initial system (charged particles in a magnetic field), so that in the 
overall system the equilibrium at rest (x l = 0) is unstable. Any deviation from this 
state, however small, is sufficient to lead to a rapidly increasing displacement along 
the peprendicular direction (Landau and Lifshitz 1976). Further analysis shows that, 
in any of the above propagation cases, the overall result would be an exponential 
increase of the perpendicular energy of the particles involved in this interaction, as a 
function of the affine parameter and, hence, of their proper-time as well. 

However, in this case, except from the primary resonance, there are also many 
secondary resonances, each one associated to a non-zero resonant width (Landau 
and Lifshitz 1976). This width is given in terms of the dimensionless amplitude 
of the gravitational wave, a, and the dimensionless ratio, (5 ~ q/uj 2 , of the energy 
density of the interstellar plasma over the energy density of the magnetic field (since 
l) 2 ~ E magn ). We may also calculate the total measure of the resonant widths, which is 
related to the total probability for our dynamical system to be trapped in a resonance. 
This is a positive quantity, directly proportional to the dimensionless amplitude of 
the gravitational wave. 

Numerical results, in the quasi-parallel case, indicate that an overall increase in the 
perpendicular energy, Ii, is also observed in any of the secondary resonances. In fact, 
in some cases this is true even when the value of Jo, in terms of which we express each 
resonant condition, lies outside of the corresponding resonant width. This indicates 
that the efficiency of the resonant interaction mechanism between charged particles 
and a gravitational wave is considerably larger than the one obtained in Papers I and 
II for the following reason. 

As a gravitational wave propagates into a non-empty space, it encounters inter- 
stellar clouds with different characteristic temperatures, densities and magnetic fields. 
To each cloud corresponds a specific set of (3, v and One expects that, for some of 
these sets, one of the many resonant conditions will be satisfied, even approximately, 
leading to a parametric resonant acceleration of the charged particles. 

We may give an estimate of the gravitational energy flow lost, due to paramet- 
ric resonant interaction of the gravitational wave with the interstellar matter, with 
respect to the flow expected at the Earth in the absence of resonant damping mecha- 
nisms. The numerical applications indicate that, for a binary source at a distance of 
30 Kpc, the total loss may be up to 10% of the flow theoretically expected. Of course 
the calculations are somewhat rough and the results refer to the idealized situation 
in which the wave propagates always parallel to the direction of the magnetic field. 
However, they indicate that the resonant interaction of a gravitational wave with the 
interstellar medium could lead to a reevaluation of the today expected values of the 
wave's amplitude in the neighbourhood of the Earth. 

Finally, the analytic results obtained in the present paper generalize the corre- 
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sponding ones of Papers I and II, since no averaging technique has been used. We 
note that the main result of the averaging procedure was to smooth out an additional 
secularily increasing term and the trigonometric modulation, which arise in the exact 
expression for the perpendicular energy h(X). 
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